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ABSTRACT 

We present an investigation of Mgll absorbers characterized as single-cloud 
"weak systems" (defined by W r (2796) < 0.3 A) at z ~ 1. We measured col- 
umn densities and Doppler parameters for Mgll and Fell in 15 systems found in 
HIRES/Keck spectra at 6.6 km s -1 . Using these quantities and Civ, Lya and 
Lyman limit absorption observed with the Faint Object Spectrograph on the 
Hubble Space Telescope (resolution ~ 230 km s" 1 ) we applied photoionization 
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models to each system to constrain metallicities, densities, ionization conditions, 
and sizes. 

We find that: 

1) Single-cloud weak systems are optically thin in neutral hydrogen and may 
have their origins in a population of objects distinct from the optically thick 
strong Mgll absorbers, which are associated with bright galaxies. 

2) Weak systems account for somewhere between 25% to 100% of the z < 1 
Lya forest clouds in the range 10 15 8 < iV(Hl) < 10 16 8 cm" 2 . 

3) At least seven of 15 systems have two or more ionization phases of gas 
(multiphase medium). The first is the low-ionization, kinematically simple Mgll 
phase and the second is a high-ionization Civ phase which is usually either 
kinematically broadened or composed of multiple clouds spread over several tens 
of km s -1 . This higher ionization phase gives rise to the majority of the Lya 
absorption strength (equivalent width), though it often accounts for a minor 
fraction of a system's JV(Hl). 

4) We identify a subset of weak Mgll absorber, those with 
logiV(Fell)/iV(Mgll) > -0.3, which we term "iron-rich". Though there 
are only three of these objects in our sample, their properties are the best 
constrained because of their relatively strong Fell detections and the sensi- 
tivity of the iV(Feii)/iV(Mgii) ratio. These clouds are not a-group enhanced 
and are constrained to have sizes of ~ 10 pc. At that size, to produce the 
observed redshift path density, they would need to outnumber L* galaxies by 
approximately six orders of magnitude. The clouds with undetected iron do not 
have well-constrained sizes; we cannot infer whether they are enhanced in their 
a-process elements. 

We discuss these results and the implications that the weak Mgll systems with 
detected iron absorption require enrichment from Type la supernovae. Further, 
we address how star clusters or supernova remnants in dwarf galaxies might 
give rise to absorbers with the inferred properties. This would imply far larger 
numbers of such objects than are presently known, even locally. We compare the 
weak systems to the weak kinematic subsystems in strong Mgll absorbers and to 
Galactic high velocity clouds. Though weak systems could be high velocity clouds 
in small galaxy groups, their neutral hydrogen column densities are insufficient 
for them to be direct analogues of the Galactic high velocity clouds. 

Subject headings: quasars — absorption lines; galaxies — evolution; galaxies — 
halos; galaxies — intergalactic medium; galaxies — dwarf 
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1. Introduction 

Absorption lines from intervening galaxies in quasar (QSO) spectra provide a wealth of 
information about the physical conditions of gas in these galaxies. Since QSO absorption 
line spectroscopy offers unmatched sensitivity to high redshifts and low column densities, the 
technique can be used to follow the gas phase in galaxies over cosmic time, from primordial 
galaxies to the local universe. 

The resonant Mgll AA2796, 2803 doublet has been used extensively to find low ionization 
QSO absorption systems (e.g. Lanzetta, Turnshek, & Wolfe (1987) and Steidel & Sargent 
(1992)). At z ~ 1, Mgll absorbers with rest frame equivalent widths, 1^(2796), greater 
than 0.3 A are optically thick in neutral hydrogen; they are observed to give rise to Lyman 
limit breaks (Churchill et al. 2000a). Furthermore, these "strong" systems almost always 
arise within 40/i _1 kpc of bright (L > 0.05 L*), normal galaxies (Bergeron & Boisse 1991; 
Steidel, Dickinson, & Persson 1994; Steidel 1995). The gas kinematics are consistent with 
material in the disks and extended halos of galaxies (Lanzetta & Bowen 1992; Petitjean & 
Bergeron 1990; Churchill et al. 1996; Charlton & Churchill 1998; Churchill & Vogt 2001). 

The discovery of weak systems (those with W r (279Q) < 0.3 A) in high resolution 
HIRES/Keck spectra (Churchill et al. (1999); hereafter Paper I) necessitates a revision in the 
standard picture. Weak systems comprise ~ 65% of Mgll selected absorbers by number, yet 
only 4 out of 19 whose fields have been imaged have a > 0.05 L* galaxy candidates within 
~ 50/i -1 kpc 5 (C. Steidel private communication). The large redshift path number density 
of weak systems relative to that of Lyman limit systems (LLSs) statistically indicates that, 
unlike strong absorbers, the majority of weak systems arise in optically thin neutral hydrogen 
(sub-Lyman limit) environments (see Figure 5 of Paper I). This has been observationally 
confirmed by Churchill et al. (2000a). Furthermore, the majority of weak systems are single 
clouds, often unresolved in HIRES/Keck spectra (resolution 6.6 km s" 1 ), in striking contrast 
to the complex kinematics of strong Mgll absorbers (Petitjean & Bergeron 1990; Churchill 
& Vogt 2001). This evidence suggests that a substantial fraction of the weak Mgll systems 
selects a population of objects distinct from the bright, normal galaxies that are selected by 
the strong Mgll absorbers. 

This begs the question: what are these single-cloud objects that outnumber Mgll- 
absorbing galaxies, yet often have no obvious luminous counterparts? A first strategy for 
addressing this question, adopted in this paper, is to constrain the column densities, metal- 
licities, ionization conditions, and sizes of single-cloud weak systems. 



5 More precisely, fields were searched within 10" of the quasar, which corresponds to ~ 50h 1 kpc at z ~ 1 
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Setting useful constraints on these physical conditions requires not just low ionization 
species, but also neutral hydrogen and medium to high ionization species. At z ~ 1, the 
Lyman series and the Cn, Civ, and Siiv transitions fall in the near ultraviolet (UV). These 
transitions were observed with the low-resolution Faint Object Spectrograph onboard the 
Hubble Space Telescope (FOS/ HST). Churchill & Charlton (1999) have demonstrated that 
photoionization modeling using both high-resolution optical spectra and low-resolution UV 
spectra can yield meaningful constraints on the physical conditions of Mgll absorbers. We 
have adapted their approach to study weak systems. 

In § 2 we describe the optical and UV spectra used in this study, and discuss sample 
selection to motivate our focus on single-cloud weak Mgll systems. In § 2.5, we describe 
our methods for Cloudy (Ferland 1996) photoionization modeling to obtain constraints on 
cloud metallicities, ionization conditions, and sizes. The resulting constraints for the fifteen 
single-cloud weak Mgll systems are presented in § 3 and summarized in § 4. In § 5 we 
consider what types of astrophysical environments could be consistent with the properties 
of the weak single-cloud Mgll absorbers. In § 6 we speculate about the evolution of weak 
Mgll systems and suggest future investigations. 



2. Data and Sample 

2.1. The HIRES Spectra 

Weak systems were charted in a survey of 26 QSO spectra (Paper I) obtained with 
the HIRES spectrograph (Vogt et al. 1994) on the Keck I telescope. A total of 30 systems 
were found, with 22 of them being new discoveries. The survey covered the redshift interval 
0.4 < z < 1.4, and was unbiased for W r (2796) < 0.3 A. The spectral resolution was R = 
45,000 (FWHM = 6.6 km s _1 ) with a typical signal-to-noise ratio of S/N ~ 30 per three- 
pixel resolution element. The survey was 80% complete for a 5<r equivalent width detection 
threshold of W r (2796) = 0.02 A. We restrict our study to those systems with this limiting 
equivalent width in the continuum at the A2796 transition. Simulations reveal that Voigt 
profile fits to HIRES data become increasing less certain below this limit (Churchill 1997; 
Churchill & Vogt 2001). This equivalent width cutoff removes from this sample three systems 
from Paper I: S9, Sll, and S22. 

Data reduction, line identification, and Mgll doublet identification have been described 
in Paper I. In Figure 1 we display all the detected transitions (Mgll, Fell, and Mgi) associated 
with the 16 weak single-cloud Mgll absorbers found in regions of the HIRES spectra that 
satisfy our equivalent width selection criterion. 
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2.2. Physically Motivated Sample 

The adopted equivalent width demarcation at 14^(2796) = 0.3 A between "weak" and 
"strong" Mgll absorbers is an artifact of observational sensitivity (e.g. Steidel & Sargent 
(1992)). However, there are at least two physical conditions that set weak systems apart 
from strong systems. First, almost all weak systems are optically thin to neutral hydrogen. 
This was shown statistically in Paper I and observationally in Churchill et al. (2000a). By 
contrast, strong Mgll absorbers are Lyman limit systems, and thus by definition are optically 
thick in neutral hydrogen. 

Second, the weak systems are often single clouds with very small velocity widths (un- 
resolved in the HIRES spectra). For the sample of 30 weak systems presented in Paper I, 
the number of clouds per weak system ranges from 1 to 7. Figure 2 shows the frequency 
distribution of the number of clouds per system for both the strong systems and the weak 
systems. The number of clouds per strong absorption system follows a Poissonian distribu- 
tion with a median of 7 clouds (see Table 7 of Churchill & Vogt (2001)). In contrast, the 
distribution for weak systems is non-Poissonian, with a spike at one cloud per absorber, i.e. 
there are many single-cloud weak systems. This suggests that these Mgll absorbers repre- 
sent a distinct population of objects. To produce many more single-cloud absorbers than 
multiple-cloud absorbers, weak systems should have a small covering factor or a preferred 
large-scale geometry, such that a line of sight is unlikely to intersect multiple clouds. 

In this paper we study the single-cloud systems. This selection criterion, together with 
the equivalent width cutoff described in § 2.1, yield 16 single-cloud systems from the sample 
of 30 systems of Paper I. 

2.3. The FOS Spectra 

For 13 of the 16 single-cloud systems, archival FOS/ HST spectra were available. The 
resolution of these spectra was ~ 230 km s _1 fwhm, with four diodes per resolution element. 
For PKS 0454 + 039 and PKS 0823-223, which contain three single-cloud weak systems, the 
spectra were retrieved from the archive and reduced in collaboration with S. Kirhakos, B. 
Jannuzi, and D. Schneider (Churchill et al. 2000a), using the techniques and software of the 
HST QSO Absorption Line Key Project (Schneider et al. 1993). The remaining six QSOs 
were observed and published by the Key Project (Bahcall et al. 1993, 1996; Jannuzi et al. 
1998) and have kindly been made available for this work. Further details of the FOS/ HST 
observations and the analysis were presented in Churchill et al. (2000a). In general, for the 
single-cloud weak systems discussed here, the useful FOS transitions were Civ, Lya, and the 
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Lyman limit, since they are strong features commonly covered in the archive data. Table 1 
lists the equivalent widths and limits for Civ and Lya, taken from Churchill et al. (2000a) 6 . 

In Figure 3, spectra covering Mgll, Fell, Civ, Lya and the Lyman limit are displayed for 
each system to show, in a single view, the full sample including blends, non-detections, and 
non-coverage. Additional transitions observed in the FOS data but not plotted in Figure 3 
can be found in Churchill et al. (2000a). 

2.4. Column Densities 

For the HIRES spectra, the column densities and Doppler parameters were obtained 
using Minfit, a x 2 minimization Voigt profile fitter (Churchill 1997; Churchill & Vogt 2001; 
Churchill, Vogt, & Charlton 2001). The column densities and Doppler parameters of the fits 
to Mgll and Fell are listed in Table 1, as are 3a upper limits on the Fell column density 
for cases where no Fell transitions were detection. The latter were obtained from the 3a 
equivalent width limit for Fell A2600, assuming a Doppler parameter equal to that obtained 
for Mgll in the same system 7 . 

The column density ratio, iV(FeIl)/iV(MgIl), is critically important to the photion- 
ization modeling; as discussed in § 2.5.1, this ratio constrains the ionization parameter. 
Therefore, it is important to appreciate the systematic errors when iV(Mgll) and iV(Fell) 
are comparable in strength. Simulations reveal that a 5a equivalent width detection limit 
of 0.02 A in HIRES spectra gives a 99% completeness limit of logiV(Mgll) = 11.9 cm~ 2 
and logTV(Fell) = 12.4 cm" 2 (Churchill 1997; Churchill, Vogt, & Charlton 2001). Above 
these column densities Minfit accurately models the column densities. For logiV(Fell) < 
12.4 cm -2 , Minfit statistically overestimates logiV(Fell) by up to ~ 0.3 dex (Churchill 
(1997), Figure 4.8) due to a bias toward "false detections" in "favorable" noise patterns. 
Mgll Doppler parameters were well recovered in the simulations, with an RMS scatter of 
~ 1 km s _1 for all column densities above the sensitivity cutoff. Fell Doppler parameters 
were poorly recovered when the column density was near or below the sensitivity cutoff 
(Churchill 1997). 

For an Fell detection at the 3a level in the lowest S/N spectrum in our sample, the 



6 Weak systems in Paper I were numbered in increasing redshift order; we have adopted those system 
numbers. 

7 Limits were insensitive to the assumed Doppler parameter because Fen was on the linear part of the 
curve of growth. 
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detection limit (99% completeness) for the ratio 7V(Fen)/7V(Mgll) would be 

log[7V(Fell)/7V(Mgll)] = 12.2 - logiV(Mgll), (1) 

Therefore, with increasing iV(Mgll) the upper limit on iV(Fell)/iV(Mgll) decreases, i.e. be- 
comes more stringent. We show this relation in Figure 4 as a diagonal line. The value of 
logiV(Fell)/iV(Mgll) for each of the single-cloud weak systems is plotted, with detections 
as solid circles and upper limits as downward-pointing arrows. Note that since the plotted 
iV(Fell)/iV(Mgll) detection limit line was computed from the noisiest spectrum in our sam- 
ple, the data points are often lower than the limit, and thus are more constraining of the 
Fell to Mgll ratio than the line indicates. 

Since there is an apparent gap in the distribution of logiV(Fell)/iV(Mgll), we somewhat 
arbitrarily define systems with detectable Fell and logiV(Fell)/iV(Mgll) > —0.3 (systems 
S7, S13, and S18) as "iron-rich" weak systems. 

Since all three so-called "iron-rich" systems lie above the worst-case sensitivity line, 
and all other systems fall below it, there may be a selection effect at work - are the high- 
est iV(Mgll) systems identified as "iron-rich" systems simply because their Fell is easier 
to detect? We argue that this is not the case. To illustrate, we consider the five sys- 
tems with iV(Mgll) between that of S18 and S7 (two iron-rich systems.) If these five sys- 
tems had logiV(Fell)/iV(Mgll) ~ 0, as is true for the iron-rich systems, then Figure 4 
makes clear that they should have detected Fell. Only one of the five systems does (S28), 
and since it has logiV(Fell)/iV(Mgll) ~ —0.6, it is not deemed iron-rich. The limits on 
logiV(Fell)/iV(Mgll) for the other four systems are also well below 0, by more than half 
a dex. Thus, it seems that while we have insufficient information to address whether the 
distribution of iV(Fell)/iV(Mgll) is continuous or bimodal, the iron-rich systems do appear 
to have significantly higher iV(Fell) to iV(Mgll) ratios than do the other systems. As § 2.5.1 
will show, this difference in column density ratios indicates variations in ionization and/or 
abundance pattern between the "iron-rich" systems and the other systems. 



2.5. Photoionization Modeling 

We assume that the single-cloud, weak systems are in photoionization equilibrium and 
constrain their properties with the photoionization code Cloudy, version 90.4 (Ferland 1996). 
The clouds are modeled as constant density, plane-parallel slabs and are matched to the Mgll 
and Fell column densities measured from the HIRES spectra. Using each model's output 
column densities for ions of interest, we synthesized FOS/ HST spectra, which we directly 
compared to the observed Civ and Lja profiles and to the spectral region covering the 
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Lyman limit break. This procedure constrains both the ionization conditions and gas-phase 
metallicities. Further discussion of the modeling technique was presented by Churchill & 
Charlton (1999) and by Charlton et al. (2000). 

We begin each analysis with the assumption that the Mgll, Fell, Lya, and Civ arise 
in a single isothermal structure, described by a single metallicity and ionization parameter. 
The data often show that this assumption is violated, in which case we model two phases, 
each having its own metallicity, temperature, and ionization parameter. Further details are 
presented below. 

Clouds were assumed to be ionized by a Haardt & Madau background spectrum (Haardt 
& Madau 1996). For cloud redshifts below z = 0.75, we used a spectrum shape and nor- 
malization at z = 0.5 and for redshifts above z = 0.75 we used the shape and normalization 
at z — 1.0. For all models, we assumed a solar abundance pattern. In § 2.5.5, we discuss 
possible abundance pattern variations and the effect of alternative spectral shapes. 

2.5.1. Applying Constraints to the Mgll Cloud 

As Figure 5 illustrates, when the assumption of photoionization equilibrium holds, the 
iV(Fell)/iV(Mgll) and iV(Clv)/iV(Mgll) ratios are uniquely determined by the ionization 
parameter, which is defined as U — n T /riH, where n 7 and % are the number density of 
photons capable of ionizing hydrogen and the total hydrogen number density, respectively. 
Therefore, \ogU = logn 7 — lognn, where n 7 is set by the background spectrum. For the 
Haardt & Madau (1996) spectrum, logn 7 = —5.6 at z — 0.5 and logn 7 = —5.2 at z — 1.0. 

When iV(Mgii) and iV(Feii) were both measured, they were used to constrain the 
optimized mode of Cloudy at a set metallicity to yield U and iV(Hi) by varying both iV(Hi) 
and n H - Models were run for a range of metallicities, —2.5 < Z < 0, in increments of 0.5 
dex 8 . 

When only an upper limit on TV(Fell) was measured, we created a grid of Cloudy models 
over log?7 (from —2 to —5, in 0.1 dex intervals) and Z (from to —2.5, in 0.5 dex intervals.) 
JV(Mgll) provided the constraint, and iV(Hl) was the parameter for which Cloudy solved. 
We rejected ionization parameters whose models, over the whole metallicity range, yielded 
TV(Fell) that was greater than the 3<r limits listed in Table 1. This set lower limits on the 
ionization parameter. 



We use the notation Z = \og(Z/Z Q ). 
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A model was judged to have failed or be inapplicable for any of the following three 
reasons: failure to converge upon an ionization parameter (usually when jV(Feii)/iV(Mgii) ~ 
1); cloud size exceeded the Jeans' length and was thus unstable to collapse; cloud size 
exceeded 50 kpc. 

2.5.2. Metallicities and Multiple Ionization Phases 

Whereas the ionization parameter is constrained by the Mgll and Fell HIRES data, 
the metallicity and JV(Hi) are constrained by the Lyct profiles and presence or absence 
of the Lyman limit break in the FOS spectra. In the regime where a cloud is optically 
thin at the Lyman limit, to create more W T (hya), given a measured iV(Mgii), requires 
lowering the metallicity, which increases iV(Hl) the total hydrogen column density Nn, where 
N n = iV(Hl) + iV(Hn), and the cloud size. 

Because the low-resolution FOS/ HST profiles are largely dominated by the instrumental 
spread function, their column densities and Doppler parameters cannot be directly measured. 
In order to use the FOS spectra as a constraint, we created synthetic FOS spectra (infinite 
sampling and signal-to-noise ratio) using the kinetic temperature and column densities 
output by Cloudy. 

First, we assumed a priori that all detected absorption arises in one phase of gas, the 
same phase that gives rise to Mgll absorption. Using the temperature of the Cloudy model 
and the measured 6(Mgii), we solved for the turbulent component, b tur b, and computed the 
Doppler parameter for other ions by the relation 

b 2 ian = 2kT/m ion + b 2 turb . (2) 

This inferred Doppler parameter and the column density are used to generate a Voigt profile, 
which is then convolved with the FOS instrumental spread function to produce a synthesized 
spectral feature. 

To constrain each cloud's metallicity, we visually compared the synthesized and FOS 
Lja profiles for each modeled metallicity. In Figure 6, we illustrate the metallicity fitting 
procedure. Metallicities of Z = 0, —1, —2, and —2.5 are plotted; clearly, Z = and Z — — 1 
do not fit the Lyct profile nor are their equivalent widths consistent with the measured 
value, and Z = —2.5 slightly overpredicts the Lya absorption. In this example, assuming a 
single phase produces the Lya absorption, the metallicity is constrained to be in the range 
-2.5 < Z < -1. 

For the systems for which the Lyman limit break is known to be absent, an additional 
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metallicity constraint may be imposed. From the Lya curve of growth, shown in Figure 7, 
for a given iV(Hi) we can find the 6(Hl) required to observe a given W(Lja). The lack 
of a Lyman limit break provides a direct constraint, logiV(Hl) < 16.8 cm -2 , which then 
implies a lower limit on 6(Hi). If this 6(Hi) is larger than that implied by thermal scaling 
of 6(Mgii) for the cloud, then there is evidence for a second phase of gas. Much of the Lya 
equivalent width must arise in this second phase, whose larger Doppler parameter produces 
the observed W r (Lja) without exceeding the limit on iV(Hi) (imposed by the Lyman limit). 
With such a second phase present, the iV(Hl) would be smaller in the Mgll cloud phase, 
and the metallicity "constraint" should be taken as a lower limit. 

The strength of the Lya profile depends only weakly on the ionization parameter for 
most of the range under consideration. Accordingly, when the ionization parameter was 
poorly constrained, we constrained the metallicity at the low and high ends of the permitted 
U range, so that the quoted metallicity constraint holds for the range of possible ionization 
parameters. 

2.5.3. Multiple Ionization Phases Required by Civ 

Often, the single-phase Mgll cloud model could not account for the observed Civ 
strength, even for the highest possible ionization parameter. If Fell is detected, it imposes 
a strict upper limit on the ionization parameter, and therefore on iV(Clv). If Fell is not 
detected, a high ionization parameter and therefore a large iV(Civ) is possible. However, 
with a small fe(Clv), scaled from the measured fo(Mgll), it may still not be possible for the 
model to reproduce the observed W r (Cw). Single-phase models were rejected when the 
equivalent widths of the synthesized profiles fell at least 3a below that measured in the FOS 
data; a second phase with a larger 6(Civ) was then required. This second phase should be 
sufficiently ionized that it does not produce a detectable broad Mgll component. However, 
its other properties (eg. its Doppler parameter and metallicity) are not well constrained by 
the low resolution FOS spectrum. 

The Civ AA1548, 1550 doublet ratio can also indicate that Civ absorption occurs in a 
broader phase. If all transitions arose in the narrow Mgll phase, the Civ doublet would be 
saturated and unresolved, with a doublet ratio of ~ 1. By contrast, a broad phase would 
produce less saturated lines, and a Civ doublet ratio closer to the natural value of 2. This 
argument is not model-dependent but relies only on the spectral resolution of FOS/ HST. 
Unfortunately, the errors in the doublet ratio are large; S17 is the only system for which the 
doublet ratio can be used to infer that the Civ arises in a broader phase than the Mgll. 
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2.5.4- Inferred Sizes and Masses 



The derived cloud size (plane parallel thickness) is simply d = Nn/nu, where Ah is the 
total hydrogen column density, iV(Hi) + iV(Hll). Assuming a spherical geometry, we can 
estimate the cloud masses within a geometric factor as 



Because % = n~,/U, this equation can also be written in terms of the ionization parameter 
where n 7 is slightly redshift dependent. 



Much of what can be inferred and/or deduced about the physical and cosmological 
properties of single-cloud weak systems is directly related to the size and mass measurements 
from the models. Here, we examine the sensitivity of these quantities to our modeling 
assumptions. 

First, we have assumed photoionization equilibrium. Photoionization models have been 
shown to underestimate the sizes of some Lya forest clouds by several orders of magnitude 
(Haehnelt, Rauch, & Steinmetz 1996). In these cases, the gas is not in full thermal equi- 
librium because of additional heating from shocks, and errors in temperature propagate to 
large errors in derived size. Such non-equilibrium is found to occur for lognn < —4 cm~ 3 , 
corresponding to the high ionization conditions of higher redshift Lja forest clouds, where 
recombination timescales can rival a Hubble time. At z — 1, such a density corresponds 
to \ogU = —1.4. As shown in Figure 5, photoionized clouds with iV(Feil) ~ JV(Mgii) are 
constrained to have logU < —3.5. Thus, iron-rich weak Mgll systems have densities too 
high, or equivalently ionization parameters too low, for their sizes to be underestimated in 
this manner. 

Second, we have assumed that the gas is ionized by a Haardt & Madau (1996) extra- 
galactic UV background spectrum. Since high luminosity counterparts are apparently rarely 
associated with weak systems, the most likely stellar contribution to the spectrum would be 
that from a single star or small group of stars quite near to the cloud. However, the con- 
straints on the stellar types, number of stars and their distance from the cloud can be quite 
severe (e.g. Churchill & Le Brun (1998), Appendix B). To explore model sensitivity to the 
spectral shape, we produced a Cloudy grid, similar to that in Figure 5, for a stellar spectrum 
characteristic of T = 30, 000 K stars. Such a spectrum, though it is too soft to ionize carbon 
into Civ, still has Fell/Mgll that decreases with increasing ionization parameter. Thus a 




M, 



.0. 



(3) 



2.5.5. Robustness to Assumptions 
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low ionization parameter (and high density) is still required for the Mgll clouds, especially 
those with measured Fell. 

Third, we have assumed a solar abundance pattern with no depletion onto dust grains. 
Enhancement of a-process elements relative to iron would shift the Cloudy grid in Figure 4 
down, so that a given iV(Feii)/7V(Mgn) ratio would correspond to a lower ionization parame- 
ter and smaller cloud size. By contrast, iron enhancement would make the clouds larger than 
we infer, but such enhancement is rarely found in astrophysical environments (Edvardsson 
et al. 1993; Mc William 1997). Dust depletion would mimic a enhancement, since iron de- 
pletes more readily than magnesium by as much as 0.5 dex (e.g. Lauroesch et al. (1996) and 
Savage & Sembach (1996a)). Dust would also affect the cooling function in the models. To 
investigate model sensitivity to this, we ran Cloudy models for S7 (a high iV(Fell)/iV(Mgll) 
system) both with no dust and with varying amounts of dust scaled relative to the ISM 
level (taken from the Cloudy database). Dust had a negligible effect on cloud sizes for dust 
abundances up to ten times ISM. At higher levels, the cloud sizes decreased with increasing 
dust content. 

To summarize, the model cloud properties are robust to the modeling assumptions. In 
particular, the conclusion seems robust that clouds with iV(Fell) ~ iV(Mgll) and low iV(Hl) 
are small. If some unknown effect has led us to underestimate the cloud sizes of the iron-rich 
clouds by 2-3 orders of magnitude, then the discussion of their origins in §5.4.1 does not 
apply, and the discussion in §5.4.2 for clouds with lower iV(Fell)/iV(Mgll) would be more 
suitable. 



3. Individual System Properties 

For each single-cloud weak system, we first summarize in brackets the coverage and 
detection status of Fell, Lya, the Lyman limit, and Civ. We then explain the constraints 
these transitions impose. Lastly, in brackets we summarize the evidence or lack of evidence 
for multiple phases of gas. Systems are listed in redshift order, with system numbers adopted 
from Paper I. The constraints for all systems are summarized in Table 2. 

3.1. SI (Q1421 + 331; z abs = 0.45642) 

[Fell not covered; no FOS spectra.] Since Lya and the Lyman limit were not covered, 
the metallicity of this system cannot be constrained. The ionization parameter cannot be 
constrained because Fell and Civ were not covered. We do not include this system in 
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discussions of multiple phases and Fell statistics because of the lack of spectral coverage. 
[Cannot address multiphase.] 

3.2. S3 (Q1354 + 195; 2 abs = 0.52149) 

[iV(Feii) upper limit; W r (Lya) measured; no break at Lyman limit; W r (Civ) upper 
limit.] To produce the strong Lya detection in this cloud requires Z < —2.5, assuming that 
all the Lya absorption arises in the same phase of gas as the detected Mgll. However, the 
lack of a Lyman limit break requires logiV(Hl) < 16.8 cm -2 , which corresponds to Z > —1.5 
for this cloud. Therefore, S3 possesses a second phase of gas with a larger Doppler parameter 
than that of the Mgll phase, which can produce the observed Lya equivalent width without 
exceeding the Hi column density limit imposed by the Lyman limit. Fell and Civ limits 
do not constrain U. [Multiphase required because of Hi; Civ does not require or rule out 
multiphase.] 

3.3. S6 (Q0002 + 051; z^ s = 0.59149) 

[iV(Fell) upper limit; Lya in spectropolarimetry mode; Lyman limit not covered; W r (Cw) 
upper limit.] The Lya spectrum is not usable because it was taken in spectropolarimetry 
mode (Churchill et al. 2000a). This, combined with the lack of Lyman limit coverage, allows 
no constraints on metallicity. The Fell limit sets a fairly high lower limit on the ionization 
parameter: \ogU > —3.5. The Civ limit is too poor to restrict the ionization parameter. 
[Cannot address multiphase.] 

3.4. S7 (Q0454 + 039; z &hs = 0.64283) 

[iV(Fell) measured; H^Lya) measured; Lyman limit not covered; H^Civ) measured.] 
If we assume the Lya arises in the Mgll phase, then to match the Lya profile requires 
— 1<Z<0. [This constraint differs somewhat from Churchill & Le Brun (1998) who 
used a different version of the FOS/i/ST spectrum.] Mgll and Fell constrain the ionization 
parameter to —4.5 < logU < —3.6 for the full range of metallicities modeled, and to —4.5 < 
logU < —4.2 for the metallicity range determined above. The model thicknesses range from 
2 pc for Z = to 8 pc for Z = —1. If the cloud were enhanced in iron, both size and 
ionization parameter would increase. The large Fell/Mgll ratio requires a low ionization 
cloud which cannot produce the observed Civ or Si IV; a second phase of higher-ionization 
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gas is required. [Multiphase is not required for Hi, but is required to explain Civ and Si IV.] 

3.5. S8 (Q0823 - 223; z ahs = 0.705472) 

[iV(Feii) upper limit; Lja and Lyman limit not covered; W r (Cw) upper limit.] Neither 
Lja nor the Lyman limit were covered, providing no metallicity constraint. The Fell and 
Civ limits set lower and upper limits, respectively, on the ionization parameter: —3.6 < 
log?7 < —2.4. In order not to exceed the 3a equivalent width, logiV(Clv) < 14 cm -2 . [Hi 
cannot address multiphase; Civ is consistent with a single phase.] 

3.6. S12 (Q1634 + 706; z &hs = 0.81816) 

[iV(Fell) upper limit; W r (Lya) poor upper limit; Lyman limit not covered; W r (Clv) 
upper limit.] The available Lya coverage is pre-COSTAR and in spectropolarimetry mode, 
and therefore unusable (Churchill et al. 2000a); accordingly, the metallicity of this system 
cannot be constrained. The Fell limit requires logC/ > —4.4. The 3<x Civ limit, which sets 
log iV(Civ) < 15 cm -2 , does not restrict U. [Hi cannot address multiphase; Civ is consistent 
with a single phase.] 

3.7. S13 (Q1421 + 331; z &hs = 0.84325) 

[iV(Fell) measured; no FOS spectra.] Without FOS spectra for this quasar, the metallic- 
ity cannot be constrained. iV(Feii)/iV(Mgii) is greater than unity for this cloud: logiV(Fell) = 
13.47 ± 0.07 cm -2 and logiV(Mgii) = 13.1 ± 0.1 cm -2 . As Figures 4 and 5 demonstrate, 
this cannot occur for a solar abundance pattern and a Haardt & Madau (1996) ionizing 
background. The iron fit is particularly robust because five Fell transitions were detected 
at relatively high S/N. Unresolved saturation may be present in the Fell transitions, which 
would cause the fit to underestimate TV(Fell). If the Mgll doublet does not have unresolved 
saturation, then it could be that this cloud is iron enhanced, as depletion cannot yield this 
pattern. A model with iron enhanced by 0.5 dex predicts logU = —4.3 and cloud size from 
1-12 pc. [Without FOS spectra, multiphase cannot be addressed.] 
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3.8. S15 (Q0002 + 051; z ahs = 0.86653) 

[iV(Fell) upper limit; W r (Lja) measured; no break at Lyman limit; W^Crv) upper 
limit.] If the observed Lya absorption were produced in the Mgll cloud, this would require 
a cloud metallicity of —2.5 < Z < —2.0, which corresponds to 18 < logiV(Hl) < 20 cm -2 . 
However, the small Doppler parameter of the Mgll cloud predicts an unresolved Lya nar- 
rower than observed, and the lack of a break at the Lyman limit 9 requires logiV(Hl) < 
16.8 cm -2 and Z > — 1 for all ionization parameters. Hence, much of the Lya absorption 
arises not in the Mgll cloud, but in a separate phase of gas with a larger Doppler parameter 
(see Figure 7). The Fell limit constrains logU > —3.6 in the Mgll cloud. Models can be 
made to meet but not exceed the Civ limit; therefore it does not constrain the ionization 
parameter; logiV(Clv) < 16.5 cm~ 2 could exist in the Mgll cloud. [Multiphase is inferred 
from Hi, though Civ does not require a second phase.] 

3.9. S16 (Q1241 + 176; z &hs = 0.89549) 

[iV(Fell) upper limit; W r (Lya) measured; Lyman limit not covered; W r (Clv) upper 
limit.] If the Lya and Mgll absorption arose in the same phase, then Lya would be best 
fit by —2.5 < Z < —1 (see Figure 6). The ionization parameter is constrained as —4.8 < 
logU < —2.0 by Fell and Civ limits. The upper limit of iV(Clv) from the 3a equivalent 
width limit, logiV(Clv) < 13.4 cm -2 , can be produced in the Mgll phase and thus does not 
require a second phase. However, if the Lyman limit were covered, its metallicity constraint 
might conflict with the Lya-determined metallicity, and thus necessitate a second phase. 
[Hi cannot address multiphase; Civ is consistent with a single phase.] 

3.10. S17 (Q1634 + 706; z ahs = 0.90555) 

[iV(Fell) upper limit; W r (Lja) measured; Lyman limit not covered; H^Civ) measured.] 
If all detected absorption arose in one phase, then Lya would best be fit by — 2 < Z < — 1, 
though these fits do not match the observed wings. However, the small Doppler parameter 
of the Mgll phase predicts unresolved, saturated Civ profiles with a doublet ratio of 1, 
which is inconsistent with the observed doublet ratio of 1.4 ±0.2. So while the log?7 = —2.0 
models can reproduce the Civ A1551 equivalent width (logiV(Clv) ~ 17 cm" 2 ), they cannot 



9 The apparent break in Figure 3 results entirely from a strong Mgll absorber at z — 0.8514 (Churchill 
ct al. 2000a). 
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fit both Civ A1548 and Civ A1551 at once. A second phase of gas, with a larger effective 
Doppler parameter, is required to give rise to the less saturated Civ profile. The Fell limit 
constrains the ionization parameter in the Mgll phase: logC/ > —3.4. [Hi cannot address 
multiphase; Civ requires multiphase.] 

3.11. Sf 8 (Q0454 + 039; z ahs = 0.93f50) 

[iV(Fell) measured; W r (Lyo;) measured; no break at Lyman limit; W r (Cw) upper limit.] 
This system has measured iV(Feii) > iV(Mgll), which a solar abundance pattern and Haardt 
& Madau spectrum cannot produce. However, because only one Fell transition was observed, 
the VP fit may be systematically large due to the noise characteristics of the data (more so 
than the formal errors would indicate). Models will converge for iV(Feii) at least la less 
than measured. We quote constraints for the range of JV(Fell) reduced by la to 2a. The 
Lya profile constrains the metallicity —1<Z<0; the lack of a Lyman limit break requires 
Z > — 1. [As for S7, the constraint differs somewhat from Churchill & Le Brun (1998) 
because of the use of a different version of the FOS/ HST spectrum.] The 6(Hi) from the 
curve of growth needed to produce the observed W r (Lya) (see Figure 7) is consistent with 
the measured fe(Mgll). These models also give —4.7 < logU < —3.6 for the full metallicity 
range modeled (—2.5 < Z < 0) and —4.7 < logU < —3.8 for —1<Z<0. These models 
have sizes below 2.1 pc. Larger sizes (5 to 175 pc) would result if iron were enhanced by 
0.5 dex. In the "iron-enhanced" scenario, Lya constrains the metallicity as above, but the 
Lyman limit break is slightly less restrictive, Z > —1.2, and the ionization parameter is 
higher: —3.7 < \ogU < —3.3. [Hi does not require multiphase; Civ limit is too poor to 
address multiphase.] 

3.12. S19 (Q1206 + 456; z abs = 0.93428) 

[iV(Feii) upper limit; W r (Lya) measured; no break at Lyman limit; W r (Cw) measured.] 
If all the Lya absorption arose in the Mgll phase, this constrains the metallicity to be 
— 2 < Z < — 1. This does not conflict with the constraint that logiV(Hl) < 16.8 cur 2 , as 
required by the lack of a Lyman limit break 10 . Thus, a second phase of gas is not required 
to explain Lya and the Lyman limit. [The value of 6(Hl) ~ 20 km s _1 from the curve of 
growth in Figure 7 is consistent with fr(Mgll) = 7.5 km s -1 for the model temperature.] The 



10 The apparent break is entirely accounted for by the z = 0.9276 system along this same line of sight 
(Churchill & Charlton 1999). 
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Fell limit requires \ogU > —3.7. The Civ detection can be produced in the Mgll cloud 
if logC/ = —1.7; however, this slightly underproduces Lya and requires an unreasonably 
large (80 kpc) cloud. Thus, multiphase ionization structure is required. [Hi does not require 
multiphase; multiphase is required to explain Civ.] 

3.13. S20 (Q0002 + 051; z ahs = 0.95603) 

[W r (2600) upper limit; W r (Lyct) measured; no break at Lyman limit; W r (Cw) mea- 
sured.] If all the Lya and Mgll absorption arose in the same phase, then the Lya profile 
would require Z < —2.0, which corresponds to logiV(Hi) > 17.3 cm -2 . Contradicting this 
is the absence of a Lyman limit break, which requires Z > — 1. Thus, a second phase with 
a larger Doppler parameter is needed to create the Lya absorption (see Figure 7). The 
Fell limit constrains the Mgll phase to logU > —3.7. The Civ absorption cannot arise in 
the Mgll cloud, even for logU = —2.0, the highest ionization parameter for which Cloudy 
can generate the observed iV(Mgll). Thus, the Civ also requires a second phase with large 
Doppler parameter, which is more highly ionized than the Mgll phase. [Multiphase is re- 
quired to explain Hi and Civ.] 

3.14. S24 (Q1213 - 003; z ahs = 1.12770) 

[iV(Feii) upper limit; no FOS data.] No HST spectra were available for this quasar, so 
the metallicity of this system cannot be constrained. We doubt the veracity of the 3.5<r Fell 
detection, since two equally strong lines are detected within 120 km s -1 of the putative Fell 
A2383 line. Were the Fell detection real, it would predict —5.2 < \ogU < —3.9 over the 
metallicity range —2.5 < Z < 0. Using the Fell as a limit, we constrain logU > —4.6. [No 
FOS data, so multiphase cannot be addressed.] 

3.15. S25 (Q0958 + 551; z ahs = 1.21132) 

[iV(Feii) upper limit; W r (Lja) upper limit; Lyman limit and Civ not covered.] Lya 
is blended with A1550 from a possible Civ doublet at z — 0.7330 (Churchill et al. 2000a), 
and the Lyman limit is not covered, so metallicity cannot be determined. The Fell limit 
constrains logU > —3.5. [Multiphase cannot be addressed.] 



-18- 



3.16. S28 (Q0958 + 551; z ahs = 1.27238) 

[iV(Feii) measured; W r (Lya) measured; Lyman limit not covered; W^Crv) measured.] 
If all the Lja arose in the Mgll cloud, —2.5 < Z < —2.0 would be required. However, the 
Lja profile is unphysically shaped and therefore this constraint is somewhat untrustworthy. 
Detected Fell and Mgll require —3.7 < \ogU < —2.8, which cannot arise in the same phase 
as the strong detected Civ. Thus, a second, more highly ionized phase is required. [Hi 
cannot address multiphase; Civ requires multiphase.] 

4. Summary of Cloud Properties 

In §3, we presented constraints on the physical conditions of fifteen 11 single-cloud weak 
Mgll systems. In this section, we summarize the inferred properties of weak Mgll absorbers. 
The inferred upper and lower limits on the metallicities, ionization parameters, densities, 
and sizes of the clouds are given in Table 2. As described in point (8) below, for at least 
seven of the systems we infer that two phases of gas are required to simultaneously fit 
the Mgll, Fell, Civ, Lja, and Lyman series absorption. In Table 2, and except where 
noted below, the inferred properties are for the low-ionization Mgll phase; the properties 
of the high-ionization phase are not well constrained because the relevant transitions were 
covered only at low resolution. Figure 8 presents an overview of the metallicity constraints 
on the low ionization phase, while Figure 9 summarizes the constraints on its ionization 
parameter /density. 

The following points characterize the basic measured and inferred statistical properties 
of the sample of single-cloud weak Mgll absorbers: 

1) Single Clouds: Two-thirds of weak Mgll absorbers are single-cloud systems, as 
contrasted with strong Mgll absorbers, which are consistent with a Poissonian distribution 
with a median of seven clouds per absorber (see Figure 2). This suggests that weak Mgll 
absorbers represent a distinct population of objects, as does their lack of Lyman breaks. To 
produce predominantly single-cloud systems, weak Mgll absorbers should have a preferred 
geometry or small covering factor if they arise in extended galaxy halos or galaxy groups, 
such that a line of sight is unlikely to intersect multiple clouds. 

2) Doppler Parameters: Most clouds are unresolved at R = 6.6 km s _1 , with Doppler 



11 There are actually 16 single-cloud weak Mgll systems in the sample. Unfortunately, no transitions other 
than Mgll were covered for SI, so photoionization models for this system cannot be constrained. While we 
include this system in number-of-cloud statistics, we exclude it from other analyses. 
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parameters of 2-7 km s 1 . A few of the Mgll profiles have slightly larger Doppler parameters 
and are slightly asymmetric, suggesting blended subcomponents or bulk motions. 

3) Metallicities: As shown by Figure 8, the Mgll phases of weak absorbers have metal- 
licities at least one-tenth solar. In no case is a cloud metallicity constrained to be lower. 
In fact, in cases where a strong Lja profile would seem to require low metallicity, assuming 
a single phase, lack of a Lyman limit break requires logiV(Hl) < 16.8 cm -2 and thus high 
metallicity. In such cases much of the Lya equivalent width arises in a broader, more highly 
ionized, second phase of gas. Thus, the Mgll phases have logiV(Hl) ~ 16 cm~ 2 and have 
been substantially enriched by metals. 

4) Ionization Parameters /Densities: Weak Mgll absorbers have a wide range of ioniza- 
tion parameters, (—5 < logU < —2) if a solar ratio of a/Fe is assumed, as shown in Figure 9. 
This translates to a range of densities, (—3.5 < logn H < cm -2 ). There is a degeneracy 
between a-process enhancement and ionization parameter. 

5) Iron-Rich Systems: iV(Fell) was detected in four of the 15 systems. Detection of iron 
considerably tightens constraints on ionization conditions (and thus sizes and densities), as 
can be seen in Table 2. Three systems (S7, S13, and S18) have iV(Fell) ~ iV(Mgll). We 
term these "iron-rich" systems, and infer logU ~ —4.5 and lognn ~ —1 cm -3 . If these 
clouds were a-group enhanced relative to solar, they could not produce the observed high 
iV(Fell)/iV(Mgll) ratio, as Figure 4 illustrates. Therefore, we infer that [a/Fe] < for these 
three systems. Also, since Fe depletes more readily than Mg (Savage & Sembach 1996b; 
Lauroesch et al. 1996), these systems do not have significant dust depletion. 

6) Systems With Fell Limits: For most of the clouds without detected Fell, we infer 
log?7 > —4 and logn H < —1-5 cm~ 3 . As Figure 4 shows, six are constrained to have iV(Fell) 
at least 0.5 dex below iV(Mgll), and thus are significantly different (either more highly 
ionized and less dense, or a-group enhanced) than the "iron-rich" clouds. For five of the 
clouds (with small iV(Mgll)) the upper limits on iV(Fell)/iV(Mgll) are not restrictive. 

7) Cloud Sizes: Together, the iV H constraint (equivalent to a metallicity constraint) and 
the density constraint (equivalent to an ionization parameter constraint) allow estimation 
of the cloud thicknesses, N-^/n-R. These constraints are given in Table 2. Iron enhancement 
would increase cloud sizes. For the three iron-rich clouds (see point 5 above), low A^h (high 
metallicity) and high density (low ionization) indicate that the clouds are small, ~ 10 pc. 
For those systems with only limits on iV(Fell), densities are not sufficiently constrained to 
infer sizes. These clouds may be as small as the iron-rich clouds, or as large as would be 
feasible for a cloud with the observed fr(Mgll) of several km s _1 (perhaps several kpcs). 

8) Multiple Phases: Seven of the 15 systems require two phases of gas: the low ionization, 
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narrow Mgll phase (which is present by definition in all systems and which may or may not 
have detectable Fell); and a second, "broader" phase that gives rise to most of the Civ 
and Lja absorption. This second phase should be rather highly ionized because broad 
Mgll absorption is not seen. The high-ionization phase is required in three systems by Lja 
and the Lyman limit, and in five systems by strong Civ. For S20, it is required by both 
Ly a /Lyman limit and Civ. 

Of the four systems with detected Fell, S7 and S28 have strong Civ which indicate 
multiphase conditions; S13 has no Civ spectral coverage, and S18 has poor Civ limits. 
Only in S18 were Lja and the Lyman limit both covered, and they do not indicate a second 
phase. 

For systems with iV(Feii) upper limits, three systems (S17, S19, S20) have Civ absorp- 
tion that require two phases, three systems (S12, S15, S16) have Civ 3a upper limits more 
stringent than the detections, and the remaining six have poor limits or no Civ coverage. 
In S3, S15, and S20, Lja and the Lyman limit indicate multiphase conditions. For all these 
systems, absorption from the second phase varies considerably in strength, indicating that in 
some cases a second phase may be absent, or at least very weak. The Lya and Civ profiles 
are of insufficient resolution to significantly constrain the properties of the higher ionization 
phase. 

9) Relationships Between Properties: Column density of Mgll, Doppler parameter of 
Mgll, metallicity, ionization parameter, and presence of a second phase are not found to be 
correlated properties in this small sample. 

5. Discussion 

5.1. High metallicity Lja forest clouds 

Numerical simulations and observations have suggested that the higher column density, 
z < 1 Lja forest clouds arise within a few hundred kpcs of galaxies (Dave et al. 1999; Cen 
et al. 1998; Ortiz-Gil et al. 1999). By stacking the spectra of 15 quasars, Barlow & Tytler 
(1998) detected Civ and derived [C/H] > —1.9 for Lja forest clouds at z ~ 0.5 (the limit 
is increased to —1.3 if clustering is considered). In contrast, at z — 3, logiV(Hi) > 15 
Lja forest clouds are seen to have a lower metallicity, with [C/H] ~ —2.5 (Songaila & Cowie 
1996). These observations indicate that larger column density, sub-Lyman limit forest clouds 
were enriched from z ~ 2 to z ~ 1. 

As Paper I demonstrated, the relative redshift number densities of weak Mgll and Lyman 
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limit systems require almost all weak systems to be optically thin in neutral hydrogen with 
logiV(Hi) < 16.8 cm -2 . For five single-cloud weak Mgll systems the FOS spectra cover the 
Lyman limit; in each of these cases there is no Lyman limit break (Churchill et al. 2000b), 
corroborating the statistical number density argument. At log iV(Hi) = 16.8 cm" 2 , for 
the measured iV(Mgll) of weak systems, the model cloud metallicities are about 10% solar. 
Unless the clouds have super-solar metallicity, logiV(Hi) cannot be much more than a dex 
below 16.8 cm" 2 . Therefore, neutral column densities of weak systems approximately range 
over 15.8 < logiV(Hl) < 16.8 cm" 2 (within a few tenths of a dex). 

We integrated the column density distribution of Lya forest clouds over this range of 
iV(Hi) to estimate what fraction of this portion of the Lya forest is associated with metal 
enriched weak systems. At z < 1.5, only the Lya equivalent width distribution is measured 
directly (Weymann et al. 1998). Converting equivalent widths to column densities introduces 
considerable uncertainty in the slope m = d\ogN/d\ogN(Rl). Using the m = —1.3 slope 
determined by Weymann et al. (1998), we find dN/dz ~ 4 for Lja over the range 15.8 < 
logiV(Hl) < 16.8 cm" 2 . Alternatively, using two power laws, with slopes mi = —1.8 and 
m2 = —0.6, which intersect at log N(HI) = 16 cm" 2 (also consistent with Weymann et al. 
(1998)), we find dN/dz ~ 1 for the same iV(Hl) range. 

By comparison, single-cloud weak systems have dN/dz = 1.1 ± 0.06 (computed as in 
Paper I). Thus, for the single power-law and double power-law fits, respectively, ~ 25% and 
~ 100% of the Lya forest clouds with 15.8 < logiV(Hl) < 16.8 cm -2 should be weak Mgll 
systems. Based upon our photoionization modeling, this implies that this portion of the 
Lja forest has been significantly metal enriched, with Z > — 1. In the iron-rich systems, it 
is likely that the Fe to Mg ratio is greater than or equal to solar, which implies enrichment 
by Type la as well as Type II supernovae (Lauroesch et al. 1996; McWilliam 1997), which 
has implications that we discuss in § 5.4.1 below. 

In Paper I, we showed that ~ 7% of the Lja forest systems with W r (Lja) > 0.1 A are 
weak Mgll systems; we have now shown that weak systems comprise a substantial fraction, 
perhaps most, of the z < 1 Lja forest with logiV(Hl) ~ 16 cm -2 . 

5.2. Space Density of Weak Mgll Clouds 

Using the sizes derived from Cloudy and the observed redshift number densities, we can 
estimate the space density of weak systems in the universe. From this we then can deduce 
the relative numbers of weak systems compared to galaxies (i.e. strong Mgll absorbers), 
regardless of how they are distributed in the universe. 
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In general, the space density of absorbers is given by 




dN\ H /c (l + 2q z) 



1/2 



dz J irR*C f 1 + z 



(4) 



where i?* is the characteristic radius of the absorber and Cf is the covering fraction, such 
that nR^Cf is the effective cross section of an absorber. We consider the Mgll phase only, 
since the high-ionization phase is poorly constrained. Since n is dependent upon the square 
of the characteristic radius of the absorber, we consider the iron-rich systems separately 
because their sizes are well constrained. 12 The systems with no measured iV(Feii) have 
dN/dz = 0.75±0.05 over the same redshift range, and their sizes are relatively unconstrained. 

Using g = 1/2, z = 1, and C f = 1, we find n = 1.3 x 10 7 /i(lpc/i?,) 2 Mpc" 3 for 
the iron-rich clouds. For R* ~ 10 pc, n — 1.3 x 10 5 h Mpc -3 (the linear dependence on 
h arises because our sizes are derived independent of cosmological model). For the clouds 
without detected iron, we find n = 5Qh (lkpc/i?*) 2 Mpc -3 . A lower limit can be estimated 
by assuming that the clouds are not Jeans unstable and are less than 50 kpc in size. This 
yields n > 0.02/i Mpc" 3 . Recall that n oc Cj 1 , so a smaller covering factor would result in 
a higher space density. 

For comparison, strong absorbers have dN/dz = 0.91 ± 0.1 (Steidel & Sargent 1992), 
R* ~ 40/i" 1 kpc, and unity Cf at (z) = 0.9, yielding n = 0.04/i 3 Mpc" 3 (Steidel, Dickinson, 
& Persson 1994), which is consistent with the space density of galaxies at this redshift (Lilly 
et al. 1995). The ratio of the space density of weak Mgll absorbers to the space density 
of strong Mgll absorbers (bright galaxies) is a simple comparison of the relative numbers 
of these populations in the universe, independent of where they arise or how they cluster. 
Therefore, though we quote the result as a number of weak Mgll absorbers per ~ L* galaxy, 
this does not imply that weak Mgll absorbers are clustered around or associated with bright 
galaxies. For the cloud size estimates of the previous paragraph, the lower limits on the 
ratio of space densities of weak Mgll systems to galaxies are ~3x I0 6 h~~ 2 and 0.5/i~ 2 for 
iron-rich clouds and clouds without detected iron, respectively. Again, recall that n for the 
weak systems goes as C^R' 2 , so a smaller covering factor and smaller sizes result in higher 
ratios. 

The sample of three iron-rich systems is small, but since space density depends only 



12 It is not clear if iron-rich weak systems are a separate population of absorbers; as Figure 4 shows, we 
cannot tell whether weak systems are bimodally or continuously distributed in iV(Feii)/7V(Mgn), and thus 
in ionization condition and/or a/Fe enhancement. For the systems with lowest ^(2796), even TV(Feii) ~ 
TV(Mgn) would place Fen below our detection threshold. Accordingly, the derived space density of iron-rich 
systems should be taken as a lower limit. 
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linearly on dN/dz (equation 4), for an unbiased survey this uncertainty does not change the 
qualitative result; more important is the robustness of their inferred size (see § 2.5.5). This 
is an astonishing result: iron-enriched single-cloud weak systems outnumber bright galaxies 
by at least a million to one. 

5.3. Relationship to Galaxies and Covering Factor 

As calculated above, the redshift number densities of all single-cloud weak systems is 
~ 1, which is comparable to that of the strong systems associated with bright galaxies. 
Taken at face value, one could argue that almost all weak systems, though significantly 
more numerous than galaxies, are nonetheless associated with galaxies. In fact, strong 
Mgll absorbers commonly have weak clouds with W r (2796) ~ 0.1 A at intermediate to high 
velocities, i.e. 40 to 400 km s -1 from the absorption systematic velocity zero point (Churchill 
& Vogt 2001). We present a direct comparison of these properties in Figure 10, which shows 
that the ranges of Mgll column densities, Doppler parameters, ionization conditions and 
iV(Feii)/iV(Mgii) are similar. 

Thirteen of the single-cloud weak systems in our sample are in QSO fields that have 
been imaged in efforts to identify the galaxies hosting strong Mgll absorption. We cite a few 
examples to illustrate that there is growing evidence against the association between bright 
galaxies and single-cloud weak systems (also see Paper I). 

Of the iron-rich systems, two (S7 and S18) are seen in absorption against Q 0454 + 039, 
whose field has been well studied (Steidel et al. 1993; Le Brun et al. 1997; Churchill & Le 
Brun 1998). No candidate galaxies are seen at the weak system redshifts within ~ 10" of the 
QSO, down to ~ O.OIL^ (Le Brun et al. 1997). The line of sight to Q 0002 + 051 has five 
Mgll systems, of which three are single-cloud weak systems (Paper I). The strongest (S6), 
with W 7 ,, (2796) = 0.29 A, is at the redshift of a bright galaxy, while the two others, S15 and 
S20, are unmatched with galaxies out to 20". The line of sight to Q 1421 + 331 has four 
Mgll systems, of which two are single-cloud weak systems (SI and S13, the third iron-rich 
system). There is one bright galaxy in the field with an unconfirmed redshift; statistically, 
it is likely to be associated with one of the two strong Mgll absorbers (Steidel, Dickinson, 
& Persson 1994; Steidel 1995). 

This evidence that weak systems do not commonly select bright galaxies within ~ 
50/i _1 kpc of the absorbing gas does not rule out either dwarf galaxies (or smaller mass 
objects) or bright galaxies within 100-200 kpc of the QSO. That is, they could arise in low 
luminosity, small mass structures that are either clustered within a few hundred kpc of bright 
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galaxies or are distributed in small groups of galaxies (analagous to the Local Group). If so, 
this might lead us to consider whether the intermediate and high velocity weak Mgll clouds 
in strong systems do not reside within ~ 50/iT 1 kpc of galaxies, but instead are the same 
objects as the weak systems; that is, the weak clouds in strong systems sampled by the lines 
of sight that select bright, Mgll absorbing galaxies could in principle be interloping weak 
systems that arise throughout the group environment. However, as we will show, there are 
problems with this scenario. 

If the single-cloud weak Mgll absorbers and the kinematic outliers of strong systems 
are the same population (with the same spatial distribution and covering factor) then they 
should have similar redshift number densities. This should be true regardless of whether 
they are clustered within 100-200 kpc of galaxies or distributed throughout groups. If they 
are not the same population of objects, then their relative redshift number densities gives 
the ratio of their absorption cross sections, 

(dN/dz) w = Ojw = {CfnRl) w 
(dN/dz) ~ a Q ~ (C f nRl) ' 

where the terms are the same as described in Equation 4, and where w and o denote weak 
systems and kinematic outliers of strong systems, respectively. 

In strong Mgll absorption systems, the chance of intercepting one or more intermediate 
or high velocity, single-weak outliers is roughly 55%; in roughly half of those systems, two or 
three single-weak clouds are observed (Churchill & Vogt 2001). This translates to a redshift 
number density of (dN/dz) Q ~ 0.7, regardless of the spatial relationship between galaxies 
and the intermediate and high velocity outliers. 

L* galaxies cover only a very small fraction of a group; consequently, passing through 
two strong systems or a strong system and a weak system should be similarly improbable, 
unless weak absorbers cluster strongly around bright galaxies. If the weak Mgll absorbers 
were clustered within 100-200 kpc of strong Mgll absorbers, a candidate L* galaxy would 
be observed, within 40/i" 1 kpc, 5-10% of the time, which does not violate the observational 
constraint. However, in this scenario, since the weak absorbers are spread over a fairly large 
region, their covering factor within that 100-200 kpc radius is also fairly small, and so one 
would only expect a high velocity weak outlier in 10-20% of strong systems. This conflicts 
with the observed fraction, which is close to 100% (some systems have multiple outliers.) 

Thus, to explain the kinematic outliers of strong systems, the weak Mgll absorbers 
would need to be strongly clustered around the strong absorbers. They would be rare at 
large distances. However, they would also be rare at small distances, or else more candidate 
galaxies for weak Mgll absorbers should have been found within ~ 50/i _1 kpc. This cluster- 
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ing, only at intermediate distances, seems rather contrived. A further complication is that 
clouds within 100 kpc should merge with the nearby galaxy on timescales of < 10 9 years. 

We therefore conclude that while weak Mgll absorbers and outliers of strong Mgll 
absorbers are similar in their physical properties, they generally arise in different types of 
hosts. In the case of the outliers, the host is apparently the bright galaxy responsible for 
the strong Mgll absorption. The high velocity clouds around this galaxy could arise from 
gas tidally stripped from companions, or through energetic events which eject gas from the 
disk (Churchill et al. 1999). The processes which give rise to the weak Mgll absorbers may 
be similar, but may occur in less massive, less luminous hosts. 

5.4. Possible Environments of Weak Mgll Absorption 

If we take the inferences presented above as reasonable approximations for the true 
physical conditions of weak Mgll absorbers, then what are the implications? In what envi- 
ronments does weak Mgll absorption arise? 

5.4-1- Iron-rich clouds 

We now consider the environments of the iron-rich weak Mgll absorbers (with dN/dz ~ 
0.2). These absorbers, with comparable column densities of Fell and Mgll, are inferred to 
be high metallicity (> 0.1 Z Q ), small [logA(Hl) < 16.8 cm -2 and size ~ 10 pc] gas clouds, 
which, if spherical, would outnumber bright galaxies by a factor of ~ 10 6 . They seem not to 
be closely associated with bright galaxies, i.e. within ~ 50/i _1 kpc of > 0.05 L* galaxies. 

The small inferred gas masses (a few M ) and small velocity dispersions (b ~ 6 km s _1 ) 
of the Mgll phase in the iron-rich systems suggest objects that would not be stable over 
astronomical timescales. Either the gas is transient or the clouds are confined by outside 
gas pressure or stellar and/or dark matter. The Mgll phase would be a condensation inside 
a larger structure, which may give rise to the higher ionization, larger Doppler parameter 
phase. However, it is important to note that these two phases would not be in simple 
pressure equilibrium, as they have similar inferred temperatures (if both are photoionized) 
and different inferred densities. 

The high metalli cities (> 0.1 solar) of weak Mgll absorbers require substantial enrich- 
ment. Yet because the gas in the iron-rich systems is not a-group enhanced relative to 
solar, it cannot have been enriched solely by Type II SNe or galactic winds produced by 
multiple Type II SNe, as material processed in this way is observed to be a-group enhanced 
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by 0.5 dex (McWilliam 1997; Lauroesch et al. 1996; Tsujimoto & Shigeyama 1998). There- 
fore, weak Mgll absorbers with high Fell must have retained enriched gas from Type la 
SNe. In general, retaining high-velocity SNe ejecta requires either a deep potential well or 
"smothering" of the explosion by gas in the surrounding medium. For weak Mgll absorbers, 
the lack of Lyman limit breaks and associated bright galaxies within ~ kpc argues 

that L* galaxy potential wells are not responsible. Accordingly, we consider how Type la 
supernova ejecta might be retained by smaller potential wells. 

In the pre-dark matter era, Peebles & Dicke (1968) explored the formation and expected 
properties of Population III star clusters. We explore the updated general scenario of a dark 
matter mini-halo of 10 6 -10 8 M Q with a virial velocity of tens of km s" 1 . Such a mini-halo 
could contain a dwarf galaxy or only a star cluster within it (Rees 2000). When the first 
massive stars in the cluster exploded as Type II SNe, the resulting superbubble would have 
driven much of the surrounding gas out into the halo. This process should destroy small 
halos, which sets a lower limit on the mass of halos which survive. In sufficiently large halos, 
the superbubble gas should be slowed as it sweeps out into the halo gas. Eventually the 
shell should slow to the virial speed of the halo, and may cool and fragment or mix with the 
halo gas and disperse. As the product of Type II SNe, this gas should be a-group enhanced. 
Condensations within the superbubble remnant might give rise to detectable Mgll but not 
Fell absorption. 

After the requisite delay time (~ 1 Gyr), Type la SNe should detonate within the star 
cluster. Scaling roughly, 10 5 M in stars should produce one Type la supernova per 10 9 yrs, 
assuming that the Milky Way SNe rate scales to lower mass structures. If this Type la 
supernova gas were retained, mixed with already a-group enhanced gas, and condensed, 
such a structure might be observed as an iron-rich, high-metallicity weak Mgll absorber. In 
the absence of a large potential well, trapping the debris would require smothering, either 
at small radii within a parent star cluster, or at larger radii within the surrounding halo. 
If the Type la ejecta were trapped within a star cluster, the observed small Mgll Doppler 
parameter would represent the low virial speed of the cluster; if the ejecta were trapped 
within the halo, the small fr(Mgll) would indicate a condensation in the supernova shell. 
Burkert & Ruiz-Lapuente (1997) have also considered the effect of Type la SNe on the gas 
in dwarf spheroidals. 

To summarize, the high inferred metallicities and lack of nearby bright galaxies imply 
that weak Mgll absorption arises in metal-enriched gas inside small dark matter halos. To 
consider what type of luminous structures could exist inside the halos (dwarf galaxies or 
star clusters) and how the absorbing gas is distributed within the halo (concentrated within 
star clusters or at large in the halo), we must balance the two factors that determine the 
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absorption cross section: number of parent halos per L* galaxy and the number of absorbers 
within these halos. Physically, the latter factor is determined by the generation rate of the 
absorbing gas and the persistence of the structure. 

Equation 5 simply relates absorber sizes and number densities of two populations to 
the ratio of their redshift number densities. Because the absorption statistics of strong 
Mgll absorbers are well-established at redshift 0.3 < z < 2.2 (Steidel & Sargent 1992), 
and because L > 0.1 L* galaxies are largely responsible for this strong absorption, it is 
useful to compare them to weak Mgll absorbers. We can rearrange equation 5 and use 
R s = 40/2T 1 kpc, dN/dz s = 0.91 ±0.1 for strong systems, dN/dz w = 0.18 ±0.01 for iron-rich 
weak systems, and unity covering factor C/. Then, n w /n s is the ratio of the number of halos 
containing iron-rich weak absorbers to the number of strong absorbers. The result is that, in 
each weak absorber halo, weak Mgll absorption covers the same area as a circle with radius 
R w = 17/i~ 1 (7i w /n a )~ 1/2 kpc. 

For the Milky Way, if only the dozen known dwarf satellites contribute to iron-rich 
weak Mgll absorption, then R w — 7 kpc for h = 0.7. Obviously, a very large fraction of each 
dwarf must give rise to the absorption in this scenario. Simulations generically predict more 
dark matter halos per poor group than are observed as dwarf galaxies in the Local Group 
(Klypin et al. 1999; Moore et al. 1999). For a typical L* galaxy, simulations by Klypin et al. 
(1999) of poor groups produce about 500 dark matter halos with v C i rc > 10 km s -1 per L* 
galaxy. Using this for 77 yields R w ~ 1 kpc per small halo for h = 0.7. Even with this large 
population of satellites, a large fraction of each halo would need to give rise to weak Mgll 
absorption with high 7V(Fell)/7V(Mgll). 

Dwarf galaxies and faint dark matter mini-halos might be expected to cluster less 
strongly than brighter galaxies. If they exist in abundance in voids, then this would raise the 
number of weak absorber halos per L* galaxy, n w /n s , and decrease the effective absorption 
radius per halo, R w . 

If Population III star clusters exist inside numerous small dark matter halos, gas trapped 
within the clusters might give rise to weak Mgll absorption. The correspondence between the 
virial velocity of a globular cluster and the small Doppler parameter of weak Mgll absorbers 
is suggestive, as is the similarity between globular cluster radii and the inferred sizes of the 
iron-rich weak Mgll clouds. Two problems of this scenario are that sufficiently small halos 
should be destroyed by the initial burst of Type II SNe, and that rogue star clusters with 
concentrations similar to that of Milky Way globulars would have been detected in the Local 
Group. More diffuse clusters might remain below present detection thresholds. The expected 
concentration is unknown, as it is difficult to calculate the expected packaging of the first 
and second generation of stars to form in the Universe (Rees 2000; Abel et al. 1998, 2000). 
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St ill , unless there are more than a million mini-halo hosts of iron-rich weak Mgll ab- 
sorption for every L* galaxy in the Universe, the absorption cross section per small halo 
spans more area than the inferred 10 pc size of an iron-rich weak Mgll absorber. Assuming 
spherical geometries, this would require multiple sites which could give rise to weak Mgll 
absorption per small halo. In the picture where enriched gas is trapped within star clusters, 
this would require multiple star clusters per halo. 

Alternatively, the gas could exist not in small, isolated structures but rather in sheets 
within the halos, which would explain the small sizes inferred for the absorbers and the large 
cross section for absorption per halo. This would correspond to the general picture discussed 
where enriched gas from supernovae, which has fragmented and cooled, is trapped within a 
halo. 



5.4-2. Clouds Without Detected Fell 

The sizes of the clouds with upper limits on JV(Feii) are not constrained. They could be 
as small as the iron-rich clouds or much larger, as large as narrow, single-cloud kinematics 
permit. Like the iron-rich clouds, these clouds with smaller ratios of iV(Fell) to iV(Mgii) 
arise in Z > —1 environments. However, the clouds without detected Fell do not require 
[a/Fe] ~ 0; they could be a-group enhanced. So unlike the iron-rich clouds, clouds without 
detected Fell could be wholly externally enriched. Given that they are apparently not closely 
associated with bright galaxies, two possible origins for their high metallicities are apparent: 
external enrichment from larger structures, or trapping of local SNe ejecta. 

In the external enrichment scenario, the winds and superbubbles of large galaxies pollute 
the intragroup gas and the low mass structures in the group: the low mass galaxies, tidal 
debris, and infalling clouds. X-ray observations of poor groups indicate that such gas would 
have relatively high metallicities (> 0.1 solar) (Mulchaey 2000). The level of a-enhancement 
depends on the ability of the galaxy group to retain the early Type II supernova ejecta which 
may escape into the intergalactic medium (Davis et al. 1999; Finoguenov & Ponman 1999). 
In this case, systems with lower Fell would arise in galaxy groups, and might be thought of as 
low neutral column density (sub-Lyman limit) high velocity clouds (HVCs). These would not 
be analogous to HVCs observed locally in 21 cm emission, which have logiV(Hl) > 18 cm~ 2 
and may or may not be of extragalactic origin (Blitz et al. 1999; Charlton, Churchill, & 
Rigby 2000). Rather, they would be more like the sub-Lyman limit HVCs observed in Civ 
absorption around the Milky Way, which are likely extragalactic (Sembach et al. 1999). Such 
Civ HVCs are consistent with a highly ionized single phase, but their Sill and Cn detections 
are also consistent with what would be expected for weak Mgll absorbers. 
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Because of the degeneracy between high ionization and a-enhancement, we do not know 
whether the clouds without detected Fell are a-group enhanced. If they are not, and instead 
have solar a/Fe, then like the iron-rich clouds, they cannot have been enriched by a-group 
enhanced external gas. In this case, the argument for the high Fell clouds applies: in the 
absence of a large potential well, Type la SNe ejecta need to be smothered and trapped by 
nearby gas. In this scenario, clouds with lower Fell would be more highly ionized or a-group 
enhanced versions of the iron-rich clouds. 

If the host star cluster is coeval, like a globular cluster, then after the initial burst of 
Type II SNe, all successive SNe should be Type la. If the stars arise in something more like a 
dwarf galaxy with more continuous or stochastic star formation, then Type II SNe remnants 
as well as high-ionization pockets of Type la SNR could give rise to lower Fell systems. 
Differences in the velocity spread, state of ionization, and absorption strength of the second 
phase seen in many of the weak Mgll absorbers might reflect differing host environments. 

6. Evolutionary Histories and Further Investigations 

6.1. Weak Mgll systems at other redshifts 

How would these z ~ 1 weak Mgll absorbers appear at higher redshift? At z > 1, the 
metagalactic ionizing flux was stronger, rates of star formation and galaxy interaction were 
higher, less material had condensed into galaxies, and less time had transpired for production 
of Type la SNe than at z = 1. Thus, clouds with the same total hydrogen column density 
would have been more highly ionized, a-group enhanced, and lower in metallicity. (We 
caution that since the star formation history of these systems is unknown, it is difficult 
to predict metallicity and enhancement evolution.) Higher ionization and lower metallicity 
would make Mgll absorption weaker, perhaps below detection thresholds. Consequently, 
sub-Lyman limit weak Mgll absorbers may be rarer or non-existent at high redshift. Due 
to increased ionization and possible a-group enhancement, Fell detections for weak Mgll 
absorbers should become rarer at z > 1. 

If this population of objects would not be common at high redshift, what kinds of objects 
would be selected by weak Mgll absorption? Some weak Mgll absorbers at high redshift 
might still be below the Lyman limit, but the threshold for detecting Mgll should be pushed 
to higher Nr by metallicity and ionization effects. Because of the column density distribution 
function, this would lead to fewer absorbers per unit redshift. However, both cosmo logical 
evolution and the fact that small clouds had not yet merged into large structures would have 
the opposite effect. Regardless of the relative numbers, it is likely that at a sufficiently high 
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redshift, Lyman limit systems would be detected as weak Mgll absorbers. Thus, objects 
physically associated with weak Mgll absorbers at high redshift may be completely different 
from those at z ~ 1, and may be associated with bright galaxies. 

How would the weak Mgll absorbers at z ~ 1 have evolved to the modern epoch? Under 
the less intense metagalactic flux of z — 0, detectable Mgll absorption can arise in clouds 
with less neutral and total hydrogen than at z ~ 1. This would make weak Mgll absorbers 
more common at the present day, since they would extend further down into the Lya forest. 
However, cosmological expansion and destruction through mergers should have the opposite 
effect. Such low ionization clouds should also have higher Fell then at z ~ 1, since there 
should have been more time for Type la SNe to occur. This may not be true for absorbers 
with recently-formed dwarf galaxy hosts. 



6.2. Future Investigations 

Additional studies at z ~ 1 can test some of the inferences of this paper and further 
constrain the properties of weak Mgll absorbers. Here we briefly discuss four promising 
avenues: spectroscopy of Mgll absorbers in multiply lensed QSOs, STIS UV spectroscopy, 
searches for Civ without Mgll, and narrow-band imaging. 

1) The critical inference that high Fell clouds have sizes ~ 10 pc can be tested by 
finding weak Mgll systems in the spectra of multiply-lensed quasars. The best constraint on 
absorber sizes thus far was derived from a z = 3.6 absorption system in Q1422 + 231, which, 
due to lensing, is probed by two lines of sight separated by 13/i -1 pc (Rauch, Sargent, & 
Barlow 1999). This sub-Lyman limit system has complex absorption spread over 400 km s -1 , 
observed in both low and high ionization transitions. In the low ionization transitions (Cll 
A1334 and Sill A1260), column densities vary by a factor of up to 10 between the two lines 
of sight. In particular, the reddest component was detected in only one of the two sightlines. 
The inferred density, gas mass, and metallicity of this component are consistent with the 
inferred values for high Fell weak Mgll clouds. Based on Cloudy models, this component 
would have detectable weak Mgll absorption and iV(Fell) ~ iV(Mgll). 

Size constraints have also been determined directly for Mgll absorbers, though at larger 
spatial scales. Eight weak Mgll absorbers have been observed in the z = 3.911 QSO 
APM08279 + 5255, in a very high signal-to-noise spectrum that combined light from multiple 
images. (The two brightest images are separated by 0.35".) For 3 systems at z = 1.21, 1.81, 
and 2.04, Mgll A2796 and A2803 cannot be fit simultaneously with Voigt profiles, which 
implies that, due to partial covering of the images, the column densities are significantly 
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different along the two major lines of sight, with separations ranging from 0.5-1.5/i 1 kpc 
(Ellison et al. 1999). 

2) With low resolution UV spectra, the properties of the higher ionization phase cannot 
be well constrained. High resolution spectra with STIS/HST will soon be available for some 
of the quasars in this sample. This new data should reveal whether the Civ resolves into 
multiple components at R — 30,000 (10 km s" 1 ), and whether the high ionization phase is 
offset in velocity from the low ionization phase. If the low ionization phase is due to ejecta 
in a larger halo, we might often see a velocity offset of tens of km s" 1 . For the systems for 
which Civ was not detected with FOS, we can determine whether the high ionization phase 
is truly absent. With additional transitions, ionization conditions in the higher ionization 
phase can also be determined. 

3) Sensitive searches for Civ at z ~ 1, especially Civ with no corresponding Mgll 
absorption, would constrain the relative sizes of the Civ and Mgll phases, given the picture 
that narrow weak Mgll absorption arises in a condensation surrounded by a broader higher 
ionization phase. 

4) It would be very time-consuming to search for galaxy groups at the redshifts of 
known Mgll absorbers via wide-field imaging and spectroscopy. Narrow-band imaging is 
more feasible (Yanny & York 1992), although small, low luminosity galaxies directly in front 
of the QSO would still be missed. 

Low-redshift investigations may also be relevant. Weak Mgll absorbers at low redshift 
should be detected serendipitously in STIS/HST QSO spectra. Because of the small redshift 
path-length for detection of Mgll absorption, few detections are expected. Nevertheless, such 
detections would help to constrain the evolution of dN/dz, and at low redshift, searching 
for associated luminous structures may be more feasible. Also, deep 21 cm mapping and 
absorption studies should shed light on how high velocity and other Lyman limit clouds are 
distributed around nearby galaxies and within nearby groups, further probing the nature of 
faint, low-column density structures in the universe. 

7. Conclusion 

The basic properties of the single-cloud weak Mgll absorbers were outlined in § 4. We 
conclude the paper by summarizing our discussion of the nature of these absorbers and their 
relationship to other classes of absorbers and objects. 

1) Single-cloud weak Mgll absorbers are of high metallicity (Z > —1) and they comprise 
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a large fraction of the logiV(Hl) ~ 16 cm 2 Lya forest (see § 5.1). 

2) The physical properties of the single-cloud weak Mgll absorbers are similar to those 
of kinematic outlier clouds in strong Mgll systems. However, most weak absorbers are 
not observed within ~ 50/i _1 kpc of L* galaxies. Cross-section arguments, outlined in 
§ 5.3, indicate that the single-cloud weak Mgll absorbers and the kinematic outlier clouds 
in strong Mgll systems cannot be one and the same population of objects (viewed from 
different orientations). They can, however, have a related process of origin. 

3) Three single-cloud weak Mgll absorbers are constrained by their relatively large Fell 
column densities to have small physical sizes, < 10 pc. Their observed dN/dz, compared 
to that for strong Mgll absorbers, indicates that, if that small, they should outnumber 
L* galaxies by more than a factor of a million (see § 5.2). These iron-rich single-cloud 
weak Mgll absorbers do not correspond to any known population of object in the local 
universe. As we discuss in § 5.4.1, their Fe to Mg ratio requires in-situ enrichment by Type 
la supernovae. Their sizes and velocity dispersions suggest an origin in star clusters (the 
elusive Population III?) or in shell fragments from the supernovae. The number of iron-rich 
Mgll absorbers required is large even compared to the number of low-mass dark matter 
halos ( "failed galaxies" ) predicted by dark matter simulations. 

4) The physical properties (particularly the sizes) are not as well-constrained for the 
larger subset of single-cloud weak Mgll absorbers without detected Fell (see § 5.4.2). Unlike 
the iron-rich population these could be a-group enhanced, though their lack of association 
with bright galaxies require energetic ejection or an origin in dwarfs. The low-iron subclass 
could represent lines of sight through sub-Lyman limit regions of high velocity clouds in 
galaxy groups. Alternatively, these low-iron single-cloud weak Mgll absorbers could arise 
in fragments of Type II supernovae or in relatively high-ionization fragments from Type la 
supernovae. 

The precise nature of the objects that host single-cloud weak Mgll absorbers is not 
known. Generally, they select high metallicity pockets of material in intra-group and/or 
intergalactic space. The phase structure apparent in many of them suggests condensations 
within larger potential wells, such as dwarfs, but the large number of absorbers is surprising. 
Understanding the processes of origin of these mysterious weak Mgll absorbers is likely to 
teach us about a common, but heretofore unknown, metal-enriched class of object. 
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Table 2: Inferred Properties of Single Cloud Weak Mgn Absorbers 
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Fig. I. — The high-resolution spectra for the 16 single-cloud weak Mgll absorbers of this 
sample. Mgll, Mgl, and Fell were captured with HIRES/Keck at R =6.6 km s" 1 resolution; 
MINFIT Voigt profiles are superimposed (see § 2.4). Table 1 lists the column densities and 
Doppler parameters. 
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Fig. 2. — The distribution of the number clouds per Mgll absorption system, for a limiting 
equivalent width of 1^.(2796) = 0.02 A. The total distribution (strong and weak combined) 
is outlined, and weak clouds are shaded. Weak absorbers show a strong spike at N — 1 
cloud per system, which suggests that they are a different type of object than the strong 
absorbers. 
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Fig. 3. — The "data matrix" for the single-cloud weak Mgll systems. For each absorber the 
Mgll A2796 transition is shown in the top subpanel. In the respective lower subpanels are 
presented the spectral regions where the Fell A2600 (or A2383) transition, the Civ doublet, 
the Lya transition, and the Lyman limit break are expected. Ticks above the spectra give 
the locations where features are expected. The full velocity window of the subpanels with 
Mgll and Fell is 100 km s _1 and for the FOS data is 5000 km s _1 . "No-Cov" indicates that 
the spectral region was not observed, and "Bad-Data" indicates that signal-to-noise ratio in 
the spectral region was too low for a useful measurement. "D" indicates a clean detection 
at the 3cr or greater significance level. "L" denotes no detection, but only an upper limit on 
the equivalent width. "Bl" indicates poor constraints due to blending with other features. 
Transitions not plotted can be found in Churchill et al. (2000a). 
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Fig. 4.— The ratio of JV(Feii) to jV(Mgii) versus JV(Mgii). Clouds with detected Fell (filled 
circles) are identified by system number. Limits, depicted by open circles, are obtained 
from the 3a equivalent width limits on Fell A2600. The dotted diagonal line on the plot 
represents the 3a detection limit for Fell in a spectrum with limiting equivalent width 
W r (2600) = 0.02 A. The individual system data points are superimposed on a Cloudy grid 
for metallicity Z — — 1 with a solar abundance pattern, and using a Haardt-Madau (1996) 
spectrum at z — 1. Solid lines indicate constant log iV(Hi) and dotted lines indicate constant 
log?7. Higher metallicity would shift the grid to the right; a-enhancement would shift the 
grid down. In order that S18 does not have a significant Lyman limit break, its metallicity is 
constrained to be significantly larger than Z — — 1. High above the permitted grid of values, 
S13 may be iron-enhanced. It is clear that a solar or slightly iron-enhanced abundance 
pattern is required to produce the high iV(Feii)/iV(Mgii) ratios. 
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log U, Ionization Parameter 



Fig. 5. — The ratios iV(Feii)/iV(Mgii) and iV(Civ)/iV(Mgii) are uniquely determined func- 
tions of the ionization parameter over three dex of log U. Since ionization structure is not 
important for weak Mgll absorbers, the ratios are independent of metallicity. Note that at 
low values of log £7, the Fell/Mgll ratio flattens and thus provides less constraint. 
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Fig. 6. — Illustration of procedure to constrain cloud metallicity from Lya. Model predic- 
tions for metallicities of Z — 0, —1, —2, and —2.5 are superimposed on the Ly« profile 
of S16. Lower metallicities predict more neutral hydrogen for a given observed N(Mgll). 
Clearly, Z = and Z = —1.0 do not fit the Lya profile; thus, the metallicity is constrained 
to be -2.5 < Z < -1.5. 



-45 - 




12 13 14 15 16 17 18 19 20 

log N(HI), cm" 2 

Fig. 7. — Illustration of the constraints placed on the Doppler parameter, b, of a broad 
Lja component from the curve of growth of the Lya transition. Curves of equivalent width 
versus Hi column density are given for six different values of 6(Lyo;), ranging from 10- 
100 km s -1 . The equivalent widths of Lja are shown as horizontal lines for the five weak 
Mgll absorbers without a detected Lyman limit break. Restricting logiV(Hl) < 16.8 cm -2 
for the measured W r (2796) provides a lower limit on b(Lya) from this curve of growth. From 
these considerations, systems S3, S15, and S20 have b(Lja) > 40 km s _1 . 
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Fig. 8. — Metallicities of single-cloud weak Mgll absorbers. Metallicities apply to the Mgll 
phase only, and are constrained fairly high. Symbol type indicates whether Fell was de- 
tected, and filled symbols indicate that a second gas phase was also inferred. Both Ly« 
and the Lyman limit constrain Z. The Lyct constraint is dependent upon the assumed 
Doppler parameter, scaled from b(Mgll). By contrast, lack of a Lyman limit break re- 
quires logiV(Hl) < 16.8, regardless of Doppler parameter. When these constraints conflict, 
the cloud has two phases of gas; the Lyman limit constraint is correct and sets an upper 
limit on N(Hl), whereas the Lya constraint is in error because much of the Lya equiva- 
lent width arises in the second phase. The metallicity of the high-ionization phase cannot 
be constrained with low-resolution spectra. Metallicity only weakly depends upon U ; the 
quoted metallicity covers the range of permitted U, which was constrained by Civ limits 
and Fell limits or detections. 
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Fig. 9. — Inferred Density and Ionization Parameter. Ionization conditions were determined 
in one of two ways: 1) by detection of Fell, in which case the vertical range between the 
two datapoints is the permitted range of logU or 2) by Fell limits and less frequently by 
Civ limits, which set lower and upper limits on logf/, respectively. Ionization parameter 
and density are related by U = n^/nu, where n 7 , the number density of ionizing photons, is 
redshift-dependent. For display purposes, we have adopted logn 7 = —5.4, an intermediate 
value, to convert the derived densities into approximate log U. To obtain the actual values 
of \ogU from the plot, add 0.2 dex to \ogU values for systems S8 and below, and subtract 
0.2 dex to log?7 values for S12 and above. 
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Fig. 10. — Comparison of weak single-cloud Mgll absorbers and weak single-cloud outlying 
subfeatures (50-400 km s _1 ) of strong Mgll absorbers, for a limiting equivalent width of 
W r (2796) < 0.02 A. The top panel plots the Doppler parameter fr(Mgll) versus log N( Mgll), 
and shows little difference between the distributions for the isolated and outlying clouds. The 
lower panel shows the ratio of iV(Fell) to iV(Mgll), which is used to place constraints on 
the ionization parameter, logt7. Both the isolated and outlying weak Mgll absorbers show 
a range of Fell to Mgll ratios. For both types of systems, several objects are constrained to 
have small thicknesses, ~ 10 pc. 



